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ABSTRACT 

We revisit recent claims that there is a "cold spot" in both number counts and bright- 
ness of radio sources in the NVSS survey, with location coincident with the previously 
detected cold spot in WMAP. Such matching cold spots would be difficult if not impos- 
sible to explain in the standard ACDM cosmological model. Contrary to the claim, we 
find no significant evidence for the radio cold spot, after including systematic effects 
in NVSS, and carefully accounting for the effect of a posteriori choices when assessing 
statistical significance. 
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1 INTRODUCTION 

Cosmic microwave background (CMB) maps have been 
studied in detail during the last few years. These stud- 
ies have been motivated by the rema rkable full-sky high- 
rcsolution maps obtained by WMAP (jBennett et al.ll2003l : 
iSpcrgcl ct al . 2003), and led to a variety of interest- 
ing and unexpected findings. Notably, various anoma- 
lies have been claimed pertain ing to the alignment 
of lar gest modes in the CMB (Ide Oliveira-Costa et al.l 



20041: iHa iian k Souradcc d |2003|: [Slosar 



Scliak '2004 

Tegmark ct al.,,.2003. : Schwarz et al.ll2004 : Land fc Maguciio 



2005bl la 



Copi et al.l 2006 ), the missing power on large an- 



gular scales (|Spergel et al.|[2003l : ICopi et al.ll2007l), and the 
asymmetries in the dis t ribution of power (|Eriksen et al.l 
I2OO4I : iBernui et al]|2006l : iHaiianllioOTl ). In the future, tem- 
perature maps obtained by the Pla nck experiment, and 
large-scale polarization information (|Dvorkin et al.l |2008| ) 
may be key to determining the nature of the large-scale 
anomalies. For a review o f the anomalies and attempts to 
explain them, see iHutereij ( 20061 ) . 

Recently, a paper by iRudnick et ai] (|2007l ) attracted 
particular attention, as it claimed to have detected a 'cold 
spot" — a drop in the source density and brightness in 
the NVSS survey. This claim would be relatively unre- 
markable, if it were not for the fact that a previously re- 
ported, anomalously cold s pot in the WMAP mic r owave 
signal (IVielya et al.l |2004 (Cruz et all |2005| . I2OO6I , l2007l : 
ICavon et al.ll2005h apparently lies at roughly the same loca- 
tion. 

This claim, if verified to be true at high statistical sig- 
nificance, would represent a major result, and would be dif- 
ficult or impossible to explain in the st andard cosmo l ogica l 
mo del. One interpretatio n, proposed in llnoue &i Silkl (|2006l ) 
and lRudnick et al] (|2007l ). is the existence of a large (> 100 



Mpc) void at z ~ 1, which gives rise to an NVSS underden- 
sity directly, and gives rise to the WMAP cold spot via the 
nonlinear ISW effect. However, the probability of forming 
such a large void in ACDM cosmology is negligibly small. 



Here we reexamine these claims using our own analysis 
procedure, carefully including known systematic and statis- 
tical properties of NVSS (declination-dependent "striping" 
and galaxy-galaxy correlations; see i]3.1fl . and marginaliz- 
ing a posteriori choices when assessing statistical signifi- 
cance. We will argue that there is no statistically signifi- 
cant evidence for either a dip in NVSS number counts or 
median source fiux in the WMAP cold spot. We will see 
that it is possible to construct statistics containing a pos- 
teriori choices (e.g., the location and radius of a "sub-disc" 
of the cold spot) which might appear to support an under- 
dense region, but the statistical signifiance goes away when 
these choices are properly marginalized. Furthermore, we 
will show that by making different a posteriori choices, we 
could find evidence for an overdense region with the same 
statistical significance as an underdense region. 



The paper is organized as follows. In ij2]we describe the 
NVSS data and selection cuts that we consider. In ij3] we 
perform statistical tests using the number density distribu- 
tion of NVSS sources, and in ^we do the same for the fiux 
distribution. In ij5] we study the dependence of our results 
on the selection cuts that are applied to the NVSS catalog 
prior to the analysis. Our main result, showing significance 
of anomalous number counts or source fluxes in the WMAP 
cold spot, using several different statistics and with a range 
of possible selection cuts in the NVSS catalog, is shown in 
Tabled] We conclude in ^ 
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Figure 1. Galaxy count maps for the NVSS survey, smoothed 
to 1° radius and plotted in equatorial coordinates, with default 
data cuts as described in [[2] The full sky counts (top panel) 
show declination-dependent variations in completeness level; the 
WMAP cold spot (shown as a circle in both panels) is entirely con- 
tained within the underdense "stripe" at declination 5 < —10°. If 
we zoom in on a box at the same declination as the cold spot, then 
the WMAP cold spot does not look anomalous by eye (bottom 
panel; shown as a rectangular region in the top panel). 

2 DATA 

The NRAO VLA Sky survey (NVSS) is a 1.4 GHz contin- 
uum survey, covering 82% of the sky, with a source catalog 
containing over 1.8 x 10'' sources that is 50% complete at 
2.5 mjy. Away from the galactic plane, almost all of the 
sources are extragalactic: quasars, or AGN-powered or star- 
forming galaxies. The NVSS catalog covers a wide range of 
redshifts (the median redshift is z ~ 0.9), but dividing the 
catalog into redshift bins is not possible because per-source 
redshifts are not measured. However, in 33 we will explore 
the effect of dividing the cat alog into flux bins. 

When making Healpix jGorski et ahllioosl ) maps from 
the NVSS catalog, we mask pixels near the galactic plane 
(l&l < 10°) or the boundary of the survey (declination 
5 < —37°). Our "default cuts" will consist of this pixel mask, 
plus dropping all sources which are flagged as having com- 
plex structure, but we will also consider other choices of cuts 

in 

In Fig. [T] we show a full-sky map and a zoomed-in re- 
gion near the WMAP cold spot, with default selection cuts 
and smoothed to 1° resolution. We have shown the full-sky 
map in equatorial coo rdinates to highlight a known system- 
atic effect in NVSS (|Blake fc Walj |2002t ): the presence of 
declination-dependent variations in completeness level, most 
notably the underdense "stripe" at 5 < —10°. Because the 
WMAP cold spot is inside the stripe, modeling these varia- 
tions will play an important role in the analysis, as we will 
see in detail below. 



3 GALAXY COUNT ANALYSIS 

Let us assume that we have a serendipitous discovery of a 
underdense spot in NVSS, and that we have prior knowl- 
edge about the location and size of the WMAP cold spot. 



We assume that the NVSS underde nsity is disc- s haped with 
center P and radius r. Following ICruz et all (|2005l ). we 
take the WMAP cold spot to be a disc with center Pq 
at {l,b) = (209°,— 57°) in galactic coordinates, and radius 
ro = 5° . How are we to properly assess the signiflcance of the 
NVSS underdensity (P, r)? We consider three possibilities: 

1. P — Po,r = ro: The NVSS underdensity has the same 
center and radius as the WMAP cold spot. 

2. P = Po,r ^ ro: The NVSS underdensity has the same 
center as the WMAP cold spot but its radius is different; 
then we must assign statistical significance in a way which 
incorporates the a posteriori choice of radius. 

3. P C Po,r 7^ ro: The NVSS underdensity lies wholly 
within the WMAP cold spot but both its center and its ra- 
dius are different; then we must incorporate the a posteriori 
choice of both radius and location. 

In the following three subsections, we will consider cases 1-3 
in detail. 



3.1 Case 1: Fixed center, fixed radius 

This case (P = Po,r — ro) corresponds to the simplest 
possible question: if we count the total number of galaxies 
in the WMAP cold spot, do we get an anomalous value? We 
introduce the ratio statistic, 

Ns.i{Po,ro) 

(iVgai(Po,ro)) ^ > 

and ask whether it differs from 1.0 with statistical signifi- 
cance, where the numerator A'gai(Po, ro) denotes the number 
of galaxies in the WMAP cold spot (Po, ro) and the denom- 
inator is its expectation value. 

Two issues arise here: first, how should the denomina- 
tor (A'gai(Po, ro)) be computed? The simplest prescription 
would be to assume that the expected number density per 
unit area is equal to the full-sky NVSS average: 

(Agai (P,r-)) = 7rr^(n)fuii_,ky (2) 

where (n)fuii_sky is the mean number density per unit area 
on the full NVSS sky. This simple estimate for (Agai) is not 
satisfactory because it does not account for the underdense 
stripe (see Fig. [T|. Therefore, we also consider an improved 
prescription based on a simple stripe model. We assume that 
the expected number of galaxies in each Healpix pixel p is 
equal to the average taken over unmasked pixels p' at the 
same declination: 

(iVgal(Pr))= <^gal(p')>p'~P (3) 

pe(P,r) 

where the sum runs over pixels p in the disc (P, r), and 
(A'gai(p'))p'^p denotes the average galaxy count taken in 
pixels p' at the same declination as p. 

The second issue when studying the ratio statistic in 
Eq. IT]) is how error bars should be assigned. Here, the sim- 
plest prescription would be to assume Poisson statistics: we 
take the uncertainty in the numerator to be given by 

AiVgai(P,r) = (Agai(P,r))i/^ (4) 

This simple prescription for (AA'gai) underestimates the er- 
ror bars because it assumes that the NVSS galaxies are 
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pure shot noise, i.e. galaxy-galaxy correlations are ignoredQ 
Therefore, we also consider an improved prescription: we es- 
timate (AA'gai) directly from the data by taking the RMS 
fluctuation over alternate choices P' of ring center which lie 
at the same declination as P: 



{iVgai(P,r)>V2 



= RMSp/ 



jVgai(P',r)-(A^gai(P',r)) 
{iVgai(P',r)>V2 



(5) 



where the notation RMSp'.^p(-) denotes the RMS fluctua- 
tion taken over choices of center P' with the same declina- 
tion as P. 

With default cuts we find the following results. 

If we use full-sky averaging (Eq. ([JJ) and Poisson errors 
(Eq. Q), then the WMAP cold spot appears to be under- 
dense in NVSS sources at 3.1 a. However this is simply an 
artifact of using the full-sky average galaxy density when 
computing (A'^gai); in fact, the WMAP cold spot is con- 
tained within the underdense NVSS stripe (Fig. [TJ. If we 
improve the estimate of {Nga_i) by using isolatitude averag- 
ing (Eq. ((Sjl) then the cold spot appears overdense at 1.1a. 
This is already not statistically significant, but if we improve 
the estimate of AN^aiiP, r) using Eq. ([5}, then the overden- 
sity drops to O.Sa. We conclude that the WMAP cold spot, 
taken as a whole, is not underdense or overdense in NVSS 
sources, but modeling the NVSS "stripe" plays a crucial role 
in the analysis. 



3.2 Case 2: Fixed center, floating radius 

We next consider the possibility of an underdense region 
with the same center P = Pq as the WMAP cold spot but 
arbitrary radius r < ro- We again define a ratio statistic 



(6) 



A^gai(fo,r) 
{iVgai(Po,r)> 

and compute the expectation value {Ngi^i{Po,r)) and RMS 
deviation AA'gai (Po, J") following the discussion in the pre- 
ceding subsection. 

Results using this statistic are shown (with default cuts) 
in Fig. [2] The rightmost error bars represent our most accu- 
rate ways of computing (A'gai) and (AA'gai) (Eqs. 
we find that all points are within la of the expected level, i.e. 
no evidence for an underdensity is seen. (We note that if say, 
one value of r gave an anomalous value, then we would have 
to incorporate the a posteriori choice of r when assessing 
significance; we revisit this issue in fjS]) The left and middle 
error bars represent less accurate ways of computing (Agai) 
and (AAgai) (Eqs. Q) and are shown for comparison. 



^ In principle, this could be remedied by estimating the galaxy 
power spectrum and including it in Monte Carlo simulations, al- 
though this may be difficult in practice, due to the presence of 
long-wavelength instrumental pow er in NVSS at low flux levels 
jSmith et al.ll2007l; [Ho et al.ll2008l) which may not be accurately 
modeled by an isotropic Gaussian field. In Eq. ||5ll, we have taken 
a simpler approach by averaging over choices of disc center P' 
in the real NVSS data, rather than averaging over Monte Carlo 
simulations. 

^ As a minor point, in Eq. JSj we have normalized each value 
of A'gai by its Poisson error 

(^gal)^''^ before taking the RMS; 
this slightly improves the estimate of (AA^gai) by accounting for 
variations in (A^gai) due to striping and the pixel mask. 
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Figure 2. Galaxy counts in circles of varying radii centered 
at the WMAP cold spot location, relative to expected counts 
as in Eq. ([TJ. From loft to right, the three sets of error bars 
(all 68% C.L.) represent increasingly accurate analysis methods. 
The left errorbars assume full-sky (A'gai) and Poisson (AA'gai) 
(Eqs. |[2]l, (O). The middle and right errorbars incorporate decli- 
nation striping in (A'gai) E'nd galaxy-galaxy clustering in (AA'gai) 
respectively (Eqs. l[3}, JSj). It is seen that if the NVSS under- 
dense stripe is not modeled, then the WMAP cold spot appears 
to be anomalously underdense, but the statistical significance is 
lost when the stripe is included in the analysis. 



3.3 Case 3: Floating center, floating radius 

The resu lt of the preced i ng sub section appears to contradict 
Fig. 5 in iRudnick et"al] (|200/f ). where a statistically signifi- 
cant underdense disc of radius r'o = 1° is seen. However, this 
figure has been constructed taking the disc center Pq to be 
the point {I, b) = (207.03, -54.85) rather than the center Po 
of the WMAP cold spot which is at (Z, b) = (209, -57). 

If the choice of (Pq, Tq) had an a priori motivation, then 
we would find, using an analysis similar to i]3.1l a 2. Oct under- 
density, with our default cuts. (There are other underdense 
"subdiscs" as well, e.g. we find that the subdisc centered 
at {l,b) = (206.82, —56.4) with radius 2.5° is underdense at 
3.0(7 if a posteriori choices are ignored.) However, we see no 
a priori motivation for making such choices of (Pq, Tq), and 
we must therefore incorporate the effect of the a posteriori 
choice when calculating statistical significance. 

To assess significance fairly, we incorporate the effect of 
the choice of (Pq, 7"o) as follows. Formally, given a disc (P, r), 
define the "naive number of sigmas" of its worst underdense 
subdisc by: 



A,(P,r) 



mm 

(P',r')C{P,r) 



Agal(P^r')-(Agal(P^r')) 
AAgai (P',r') 



(7) 



where the notation min(p/ ,,/)^(p ,,)(■) means that the quan- 
tity in parentheses is minimized over discs (P', r') which are 
entirely contained in (P, r). Then the existence of the "sub- 
disc" from the preceding paragraph can be rephrased as the 
statement that Ao-(Po,ro) = —3.0, where (Po,ro) are the 
cold spot center and radius respectively. 

To assess whether this value is anomalous, we evaluate 
the same statistic ("number of sigmas of the worst under- 
dense subdisc") in an ensemble of disc-shaped regions with 
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Figure 3. NVSS flux maps, smoothed by 2° median filtering as 
described in 21 shown with default cuts in equatorial coordinates. 
As in the galaxy count case {Fig.[TJ, the full-sky map (top panel) 
shows declination-dependent variations in median flux, and the 
WMAP cold spot (shown as a circle) is not anomalous by eye 
when viewed in a "box" at the same declination (bottom panel). 



the same size and declination as the WMAP cold spot. This 
way of assessing significance accounts for both the a poste- 
riori choice of {Po,r'o) and declination-dependent striping. 
More precisely, we compute Ncr{P, ro) for an ensemble of al- 
ternate choices of disc center P with the same declination as 
the WMAP cold spot center Pq, and with radius ro = 5° . We 
find that this ensemble of values has mean (N^) = —2.9 and 
RMS error (AiV^) = 0.4, i.e. the cold spot (with iV^ = -3.0) 
is typical among discs with the same radius and declination, 
and the "subdisc" described above loses statistical signifi- 
cance after the effect of a posteriori choices is taken into 
account. 



4 FLUX ANALYSIS 

In addition to number counts, the median NVSS brightness 
was also repor t ed to be low near the WMAP cold spot in 
iRudnick et al.l (|2007h . Since brightness is roughly propor- 
tional to (source counts) x (source flux), and we have al- 
ready analyzed the source counts, our perspective is that it 
is better to separate the two issues and next ask whether 
median source fluxes in NVSS are anomalous in the WMAP 
cold spot. Considering source counts and fluxes separately, 
rather than using brightness maps, has two additional ad- 
vantages. First, it allows the analysis to proceed from the 
NVSS source catalog, thus avoiding instrumental complexi- 
ties associated with working with the NVSS images, which 
have been incorporated by the NVSS team when construct- 
ing the source catalog. Second, it avoids introducing more 
a posteriori choices which must be marginalized (e.g. in 
iRudnick et al] (|2007l ). the brightness maps are convolved 
with an 800 arcsec filter to obtain a continuous field, which 
is then median-filtered in sliding boxes with side length 3.4°; 
a "dip" is then observed at a point other than the WMAP 
cold spot center.) 

We would also like to emphasize that, if the purpose of 
this analysis is to find voids, then there is no a priori moti- 
vation for considering either brightness or source fluxes; the 



best-motivated statistics would be based on number counts 
alone. Nevertheless, in this section, we will briefly NVSS 
source fluxes in the WMAP cold spot. Our median flux anal- 
ysis will be analogous to the galaxy count case from the pre- 
ceding section; we summarize our methodology and results 
here. 

In Fig. |3l we show a flux map obtained by taking the 
median flux of all NVSS sources within 2° of each pixel. 
This median-based smoothing procedure was used because 
the NVSS flux distribution contains far outliers; if the mean 
were used instead of the median, then the map would be 
dominated by a small number of rare bright sources. Note 
that declination-dependent striping is seen in the flux map, 
as seen previously for number counts (Fig. [T]). 

Given disc center P and radius r, we deflne fi{P, r) to 
be the median flux of all NVSS galaxies contained in the 
disc (P, r), and consider ratio statistics of the form: 



(8) 



We estimate the expected median flux {fi{P,r)) directly 
from the data by averaging over alternate choices of cen- 
ter P' with the same declination as P: 



{^^{P,r)) = {^^{P',r))J 



(9) 



We also estimate the error (A/i(P, r)) from the data in an 
analogous waj|f| 

RMSp..p(iVgai(P', r))'/^ Up', r) - (^*(P', r))) 

^'^''■'^ ^ (iVg.(P,0)v- -^''^ 

We then consider three cases for the disc {P,r), as in 

1. Fixed center, fixed radius (P — Po,r = ro): In this 
case, we simply evaluate the ratio statistic in Eq. ([8]), taking 
(P, r) to be the WMAP cold spot center and radius (Po, ro). 
We find that the ratio exceeds 1.0 by 0.6a, i.e. the median 
flux of all NVSS galaxies in the WMAP cold spot is not 
anomalous. 

2. Fixed center, floating radius (P — Po,r < ro): In this 
case, we compute the ratio statistic in Eq. (|8]) for a variety 
of radii centered at the cold spot center Po (in analogy with 
Fig. [2}. We find that all values are within 1.2a of 1.0, i.e. no 
anomalous value of the median flux is found. 

3. Floating center, floating radius (P C Po,r < ro): In 
this case, we choose a subdisc (P, r) of the WMAP cold spot 
(Po,ro) which appears to give the most anomalous value of 
the ratio statistic in Eq. If we look for an anomalously 
low median flux, then we flnd that the subdisc with center 
(l,b) = (206.81,-54.75) and radius 1° appears to be low 
at 2.2a, if the a posteriori choice of (P, r) is temporarily 
ignored. To assess whether this value is really anomalous, we 
proceed in parallel with the number count analysis in t)3.3l 
we evaluate the same statistic (naive "number of sigmas" Na 
of the most anomalous subdisc" ) over an ensemble of regions 
with the same size and declination as the cold spot. We get 



^ Note that we have taken the factor (N^^i)^^'^ inside the RMS 
average, to improve the estimate of (A^) by accounting for the 
— 1/2 

scaling (A/j oc N^^^ ) expected due to variations in number den- 
sity alone. 
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Na = (2.2 ± 0.3) in this ensemble, i.e. the WMAP cold spot 
(with Nc = 2.2) is a typical member of this ensemble, and 
the low flux in the aforementioned subdisc has no statistical 
significance. 



5 ALTERNATE CHOICES OF CUTS 

We have now performed an exhaustive analysis of NVSS 
number density (i|3]) and median flux (Q in subspots of the 
WMAP cold spot, with three cases depending on whether 
the subspot location and radius are determined a priori 
or a posteriori, for a total of six analyses in all. This has 
been done using our "default cuts" from fjS] we drop NVSS 
sources flagged as having complex structure to be conserva- 
tive, but do not impose flux cuts, in order to avoid making 
an a posteriori choice of flux range. 

One pos si ble l oophole remains: in Fig. 5 in 
iRudnick et al.l (|2007h . the statistical significance ap- 
pears to be much higher if only sources with fiux 5^5 
mjy are retained. In this section, we consider the general 
question; can we get a statistically significant result if we 
use selection cuts other than our default choice? 

We divide the NVSS catalog into four fiux bins, with 
delimiting values given by {3, 5, 12} mJy; these values were 
chosen to roughly divide the catalog into quartiles. We also 
consider either dropping or retaining sources fiagged as hav- 
ing complex structure in the NVSS catalog. 

In Table [1] we summarize the results of repeating the 
six analyses considered in this paper with various sets of 
cuts. This table presents many results in compressed form 
and is organized as we now explain. 

Columns labeled "Fixed center, fixed radius" corre- 
spond to case 1 in ^Mll we simply compute the total number 
(or median fiux) of galaxies inside the WMAP cold spot, and 
report the deviation (in "sigmas") from the expected value. 
A positive sign indicates a result (either count or fiux) which 
is larger than expected; a negative sign indicates a result 
which is less than the expected value. 

Columns in Table [1] labelled "fioating center, fioating ra- 
dius" correspond to case 3 in 33H11 We evaluate the statistic 
NaiPoiTo), defined in Eq. ((Tjl to measure the "worst under- 
dense subdisc" of the WMAP cold spot (Po, '"o)- We compare 
this with the mean {N^) and RMS (AA^^) taken over an en- 
semble of regions with the same size and declination as the 
cold spot. If the quantity 



(11) 



is positive, then the worst underdense subdisc in the cold 
spot is more anomalous than the ensemble mean {Na), and 
A/linder gives the significance of the anomaly in "sigmas" . We 
define a quantity TVover by computing the significance of the 
worst overdense subdisc in an analagous way. In Table [T] we 
report either A/'under (indicated by a negative sign) or TVovcr 
(indicated by a postive sign), whichever is more significant. 
A ' — ' entry means that both TVundcr and A/'ovor are negative, 
so that both the most underdense and most overdense sub- 
disc of the WMAP cold spot are less anomalous than the 
ensemble mean. 

Finally, columns in Table [1] labelled "fixed center, fioat- 
ing radius" correspond to case 2 from ^SHU with one differ- 



ence: when reporting the significance of the most anomalous 
radius r, we incorporate the a posteriori choice of radius by 
maximizing over r as in case 3. (We omitted this step for 
simplicity in ijSHl] because with our default cuts, all choices 
of r turned out to give very typical values.) 

There are a few values in Table [1] which might be in- 
terpreted as statistically significant, e.g. the example which 
motivated this section: for the fiux range S ^ 5 mJy, there is 
a subdisc (the "fioating center, fioating radius" case) which 
has galaxy counts low at 2.7a even after accounting for the a 
posteriori choice of center and radius. However, we note that 
the statistical significance goes away when complex sources 
are dropped, or if we restrict the fiux range further. Further- 
more, one can find another value in the table which sup- 
ports an overdensity with the same statistical significance 
(the "fixed center, fixed radius" galaxy count case with fiux 
range 5 ^ 5* ^ 12 mJy). 

Analogously, for the fiux analysis, there is a 2.8(j low- 
fiux subdisc in the fiux range 3 ^ S* ^ 5 mJy, but if com- 
plex sources are dropped, we find a 2.8(t high-flux disc in 
the same fiux range. Given the large number of entries in 
Table [T] a few high-significance values such as these are ex- 
pected as statistical eventsQ Since there is no clear pattern 
to the few high-significance values, and since a high source 
density/fiux region is supported roughly as well as a low 
source/fiux region, our interpretation of Table[T]is that there 
is no evidence for either NVSS number counts or median 
source fiuxes which are atypical in the WMAP cold spot. 



6 DISCUSSION AND CONCLUSIONS 



In thi s paper, we have revisited claims from iRudnick et al.l 
(|2007l ) that there is a cold spot in the NVSS radio survey 
which is statistically s ignificant and aligned with the col d 
spot found in WMAP (|Vielva et al.ll2004l : ICruz et al.ll2005l ). 
We found no evidence for either an underdensity in NVSS 
number counts, or a region of atypical median source fiux. 

Our analysis incorporates systematic declination- 
dependent striping in NVSS, by estimating quantities such 
as (Afgai(P, r)) directly from the data via isolatitude averag- 
ing (e.g. Eq. ((Sjl). In an analogous way, we incorporate statis- 
tical clustering of galaxies in NVSS by estimating variances 
over isolatitude rings (e.g. Eq. ([5])). 

Simple, direct statistical tests, such as counting the to- 
tal number of all NVSS galaxies in the WMAP cold spot, or 
taking the median fiux of all such galaxies, do not show any 
statistically significant anomaly. (This corresponds to case 
1 in llSHll) Things only become murky when one considers 
statistics with many a posteriori choices, such an anomalous 
subdisc of the cold spot (case 3) , or a choice of selection cuts 
which maximizes the quoted significance (Table [l}. We have 
exhaustively studied many such statistical tests and argued 
that when the a posteriori choices are included in the as- 
sessment of statistical significance, there is no evidence for 
an NVSS "cold spot". 



* To make this more quantitative, given only 9 independent 
events, the likelihood that one event is anomalous at the 2. So- 
level is 2a, so getting a few such anomalous values in an anal- 
ysis such as Table [T] with many different choices of cuts is not 
surprising. 
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NVSS galaxy counts analysis NVSS median flux analysis 
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Table 1. Statistical signiflcance of anomalous NVSS number counts or median flux in the WMAP cold spot, for different NVSS flux 
ranges, and with complex sources either dropped (unparenthesized values) or retained (parenthesized) in the analysis. The six columns 
correspond to the different number count and median flux analyses that we have considered in 331"[H1 As described in JS] each entry in 
the table is either the statistical significance of a region with high source density/flux (positive sign), or low source density/flux (negative 
sign), depending on which has higher significance. An entry is marked ' — ' if there is neither a positive or negative subdisc which is more 
anomalous than the relevant ensemble mean. 



As a co ncrete example , consi der the "S ^ 5 mJy" case 
in Fig. 5 of iRudnick et all (|2007f ). which seems to show a 
~ 5a underdense region, if the errorbars are taken at face 
value. We agree with the number counts that are plotted 
in this figure, but disagree that there is statistically signifi- 
cant evidence for an underdensity. Let us illustrate this by 
following this example through the steps of this paper one 
at a time. First, if we use our most accurate prescriptions 
for the RMS error (AA'^gai) and the expected count (A'^gai) 
(Eqs. ([31, (13), then we find that the statistical significance 
drops to 3.4(t; however this ignores the effect of a posteri- 
ori choices. The center of this underdense region is not the 
WMAP cold spot center, and if we account for this choice 
(and the a posteriori choice of radius) using the method of 
i|3.3l then we find that the significance decreases to 2.7a. 

This now accounts for the a posteriori choice of sub- 
disc, and appears to give a statistically significant result, but 
we have still made an a posteriori choice of selection cuts, 
by allowing complex structure and considering only sources 
with flux S ^ 5 mJy. When viewed in the larger context of 
Table [T] it is seen that these choices maximize the quoted 
"number of sigmas" of an underdensity, and simply reflect 
the large number of possible choices of selection cuts: the 
significance goes away if the cuts are changed slightly, and 
in fact a different choice of cuts would favor an overdensity 
rather than an underdensity, with roughly the same signifi- 
cance. This last observation is perhaps the most convincing 
sign that the apparent ~ 2.7a underdensity, for a single 
choice of selection cuts, is spurious. 

For the median flux analysis O, our conclusions are 
the same: we find no evidence for atypical source fluxes in the 
WMAP cold spot, after accounting for a posteriori choices. 
There are a few choices of cuts which appear to show anoma- 
lous values (if the a posteriori choice of cuts is ignored), but 
the number of such values is consistent with statistics, and 
the cuts can be tuned to support either a region with high 
or low source density/flux with roughly the same statistical 
signiflcance (Table [T|. 

We do not see reason to give preferential treatment to a 



posteriori choi c es in the analysis and selection cuts given in 
iRudnick et al] (|2007l ). and we instead considered a range of 
analyzed quantities and selection criteria. Had there been a 
physical reason or a survey-specific r equirement for the par - 
ticular treatment of raw data used in IRudnick et aL I (|2007h . 
we would have agreed with that choice being well motivated 
or even necessary. However, since we do not see such moti- 
vation, we insist on calling all such choices a posteriori. 

Despite the null result of this paper, one should not 
be disheartened. More detailed observation of the cold spot 
region in galaxy surveys will likely improve confidence about 
the existence of any over/underdensity or lack thereof. More 
generally, new WMAP data and the eagerly expected Planck 
maps expected in a few years, combined with data from a 
variety of galaxy surveys from ground and space, will provide 
a gold mine to search for signatures of the early and late- 
universe physics in the large-scale structure and the CMB. 
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